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The inelastic scattering of neutrons from magnetic excitations in the antiferromagnetic phase of 
the langasite compound Ba3NbFe3Si20i4 is analyzed theoretically. In the calculations presented, the 
strongly coupled spin-5/2 Fe triangles are accounted for as trimerized units. The weaker interactions 
between the trimers are included within the mean-field/random-phase approximation. The theory 
is compared with linear spin- wave theory, and a model is developed which leads to good agreement 
with the pubhshed results from unpolarized and polarized neutron-scattering experiments. 

PACS numbers: 75.30.Ds, 75.10.-b, 75.25.-j 



I. INTRODUCTION 

The langasite compound Ba3NbFe3Si20i4 belongs to 
the space group P321 (space group number 150), where 
triangles of Fe'^"'" {L — S — 5/2) ions are placed in a 
lattice as indicated in Figs. [T] and [H The space group 
contains no improper symmetry elements and the other 
ions may be arranged in two different mirrored ways as 
specified by the sign of the structural chirality factor = 
±1. The Hamiltonian for the S = 5/2 spins of the Fe ions 
is assumed to be 

^ = ^ E E • + E ^t(*), (1) 

where i and j are the triangle numbers, and ^ and 
r] = 1,2,3 denote the different spins in each triangle, 
as defined in Fig. [1] The Hamiltonian for the isolated 
triangles is determined in terms of the intra-triangle in- 
teraction Ji > and an anisotropy term assumed to be 
a Dzyaloshinsky-Moriya (DM) interaction: 





(Si X S2 + S2 X S3 + S3 X Si) 



(2) 



The system is antiferromagnetically ordered below 
= 27 Kiii^ The moments are confined to lie in 
the ab plane, which property is in accordance with the 
anisotropy introduced in Eq. ^ (independent of the sign 
of D^). All triangles in a certain ab plane are identical, 
and the three spins, (S^), (S|), and (S3) in the ith trian- 
gle are making an angle of 7 = e^(27r/3) with each other 
so that (S'i) + (S|) + (Si) = 0, where = ±1 defined in 
Fig.[T]denotes the two possible orientations of the ordered 
spin triangles. The moments along a line parallel to the c 
axis rotate the angle (/) = Q-c ~ e^(27r/7) from one layer 
to the next along the c axis, where e^^ = ±1 denotes the 
chirality /helicity of the helically ordered moments. 

The magnetic excitations at low temperatures in the 
helically ordered phase have been studied experimentally 
by two independent groups: Stock et al^ determined 
the dispersion relations of the spin waves propagating 



FIG. 1: (Color online) Fe triangles are stacked on top of each 
other along the c axis in the Ba3NbFe3Si20i4 crystal. The 
line parallel to the c axis through the center of the triangles 
is a three-fold symmetry axis, and the three altitudes of one 
triangle are two-fold symmetry axes. The oxygen ions are 
placed in different positions along the paths determining the 
interplanar exchange constants J3, J4, and J5. The short- 
est super-superexchange path (and bond angles closest to tt) 
definesi the structural chirality = ±1. If this path is the 
one determining J5 as defined in the figure, the correspond- 
ing chirality would be = —1, whereas = -1-1 if this 
path is the one leading to the exchange interaction J3. The 
right part of the figure shows the ordered spins in one of the 
triangles, where the moments rotate the angle 7 = e.^(27r/3) 
for increasing site index. The upper triangle shows the case 
where = -1-1, whereas the orientation of the ordered spins 
in the lower triangle corresponds to = — 1. 



in the a*c* plane from inelastic scattering of unpolar- 
ized neutrons, whereas Loire et af^ investigated the ex- 
citations by doing inelastic scattering experiments in the 
b*c* plane with both unpolarized and polarized neutrons. 
Loire et al.'^ carried through a linear spin-wave analy- 
sis of their results from which they concluded that the 
moments were all ordered in righthanded helices, corre- 



sponding to = +1, in the 



-1 enantiopure crystal 
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FIG. 2: (Color online) The triangles within an ab plane with 
a definition of the exchange parameter J2. 



they were investigating. 

The spin system in Ba3NbFe3Si20i4 is relatively 
strongly frustrated, and the validity of the mean-field 
(MF) and the random-phase approximation (RPA), on 
which the linear spin- wave theory is based, may be ques- 
tioned. The dominant cause for frustration is the strong 
interaction between the three spins in the Fe triangles. 
In the present paper I have improved the theoretical de- 
scription of the system by treating the spin triangles as 
trimerized units placed in a mean-field due to the inter- 
actions between neighboring trimers. The properties of 
the trimers are determined accurately by a numerical di- 
agonalization of the mean-field trimer Hamiltonian, and 
the corresponding spectrum of collective trimer excita- 
tions is calculated by the use of standard RPA numerical 
techniques 1^1^ Before the numerical analysis is presented 
in Sect. Ill, it is instructive to consider the approximate 
but analytical behavior of the spin system as done in 
Sect. II. The spin-wave theory applied in this section is 
not the standard one based on the Holstein-Primakoff 
transformation, but it is straightforward to show that the 
present RPA approach and the standard spin-wave the- 
ory are equivalent to leading order (see for instance Ref. 
H). In Sect. IV, the cross section for polarized neutrons 
is calculated numerically, and the results are compared 
with the experimental results obtained by Loire et alA 
The conclusions are presented in Sect. V. 



II. LINEAR SPIN- WAVE THEORY 

A coordinate system is assigned to each spin, defined 
so that the local z axis is along the direction of the or- 
dered moment, and the y axis is along the c axis and 
common for all coordinate systems. In the spin-wave 
limit, where the spins are fully polarized, the ground 
state of the ith triangle is |000). The states are defined 
by S^^\n^n2n^) = {n^ - S)\n^n2n^), where ^ = 1,2,3. 
The interactions between the different triangles lead to a 
mean field acting only on the z components (neglecting 



any anisotropy within the ah plane) : 

^MF = ('S'lz + + -Ssz) + ^K-r . (3) 

In this section it is assumed that h ^ Ji and only terms 
to first order in J^/h and D^/h are included. The intra- 
triangle interactions imply that the ground state \g) is 
modified into 

|g) = |000)-HA(|110)-F|101) + |011)), 

35* Ji —25* 13^ sin 7 ,,, 

' - sh • 

The first-order modification of the ground state does not 
affect the ground-state energy to leading order, and it is 

E{g)^-3Sh-^S'^Jj^+3S^D^sinj. (5) 

The presence of the DM anisotropy implies a specific sign 
for 7. The product D^sin^ has to be negative, and the 
orientation of the ordered spin triangles is determined by 
the sign of according to = — sign(D^). 

In terms of the interaction parameters defined in Figs. 
[Hand [21 the exchange field is: 

h = —2S[ 2 J2 cos 7 + J4 cos (j) 

+ J3 cos((?!) + 7) + J5 cos((/) - 7)] . (6) 

(j) is the angle of rotation of S^, within the plane, when 
going from one ab plane to the next in the positive c di- 
rection. The value of (j) is found by maximizing h (mini- 
mizing the ground state energy) and is determined by 

tan0 = i?sin7, R ^ ~/'\ . (7) 

It is remarkable that a difference between J3 and Jg im- 
plies that 4> becomes non-zero with no need for an inter- 
action between next-nearest neighboring spins along the 
c axis. Experimentally;^ it is found that \(f)\ is close to 
27r/7. The relation also shows that the magnetic helic- 
ity and the orientation of the ordered spin triangles are 
intimately related, that 

Eh = sign(i?) . (8) 

The three lowest excited states are (approximately) 

\a) = (|100) + w|010) + w^lOOl)) 7^3 

\h) = (|100) -I- w^lOlO) + w|001)) /V3 (9) 

|c) ^ (|100) + |010) + |001))/\/3, 

where w = e^'^^/'^ . The corresponding excitation energies, 
K = E{a)- E{g), are 

A„ = A, = /i+|5Ji-|5D,sin7 

A^ = h+ lSJ^ - . (10) 

The collective excitations of this system may be derived 
by introducing the "boson-like" (creation) operators for 
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the ith triangle: oil\g) = |a), where a = a, 6, or c. 
Neglecting the interactions between the triangles and in- 
troducing the Fourier transforms of the operators, the 
single-trimer MF Hamiltonian may be written: 

(11) 

Within the subspace of the four lowest spin states for the 
ith triangle, we may write 




and similarly for the other spin components, except that 
the {a^,b^) and {a'^,b^) terms are multiplied, respec- 
tively, by w and w'^ in the expressions for the spin com- 
ponents of Sj, and by and w in the spin components 
of S3. To first order in Ji/h, the relative matrix elements 
are 

m%^l- X, m° = 1 + A , (13) 

in terms of the mixing parameter A defined in (|4]). Next 
step is to substitute these expressions for the spin com- 
ponents in the exchange Hamiltonian, which leads to an 
effective Hamiltonian quadratic in the excitations oper- 
ators. Introducing the Fourier transforms of the opera- 
tors, the Hamiltonian stays diagonal with respect to the 
(Cq,Cq) operators, when q is parallel to the c axis . We 
shall concentrate on this case in the following, i.e. that 
q • c = gc. In the zero temperature limit, the only non- 
zero commutator relations are [aq,aq/] ~ '^qq'' where 
a = a, b, or c. From the equations of motion, the energy 
squared of the c-mode excitations propagating along the 
c axis is then found to be 

E^M = [Ae + S{mlf {[J(Q) + 2J,] cos(gc) - 2J^}] 
X [A, + S{mlY {[J(0) - 4J2] cos(gc) + AJ^)] , (14) 

where the two exchange parameters are 

J{Q) = 2 [J4 cos -f J3 cos(0 -f 7) -I- Jg cos(0 — 7) — J2] 
J(0) = 2 [2 J2 + J3 + J4 + J5] . (15) 

J{Q) is also the parameter determining the exchange 
field, h = —SJ{Q). The remaining part of the excita- 
tion Hamiltonian leads to two w modes with mixed \a) 
and \b) state characters, and the squared energies of these 
modes are 

El,{q) = [A, -f S{m'^,f3M [A, + S{m'^y)%iq)] 

ElM = [A, + S{mlf2.A-l)\ [K + S{mlf^y{-q)\ , 

(16) 



where 

?i.j.[q) — 2J4 cos (/)Cos(gc) -I- 2J3 cos((/) + 7) cos((7c + 7) 
+ 2 Jg cos((/) — 7) cos[qc — 7) + ^2 

■Syiq) = 2J^cos{qc) + 2J^cos{qc + j) (17) 
-I- 2 Jg cos{qc — 7) — 2 J2 . 

In general 3x,y{—Q) 7^ Sx^yil), which means that the two 
modes are not symmetric around q = 0, instead we have 
tha.t E^,{q) = E^^{-q). 

Introducing the expressions for A^ and the relative 
matrix elements into Eqs. (|14l) and (fTB|) . the squared ex- 
citations energies are, to leading order in A, 

E^^iq) ^S^l- cos(gc)] [- J(Q) - 2J2] (18) 
X [3 Ji - 2D^ sin7 - J(Q) -I- 4J2 + {^(0) - 4J2} cos(gc)] 

and 

Eliq) = 52 [3y{q) - J(Q) - 2i?,sin7] (19) 
X [|Ji-3i?,sin7 + 3,(g)- J(g)] . 

Equation (ITS)) shows that the energy of the c mode E^ (q) 
vanishes linearly with q, when q — >■ 0. This Goldstone 
mode, which appears because the rotational symmetry 
in spin space around the c axis is broken in the ordered 
phase, reflects that it costs no energy to rotate the or- 
dered structure around the c axis. In the case of a c exci- 
tation propagating along the c axis, all the c components 
of the spins in a certain ab plane are moving in phase, 
whereas the w modes imply that the plane of the spins 
in a certain triangle is oscillating out of the ab plane, and 
this oscillation is affected by the DM anisotropy. Because 
of this anisotropy, the w modes show an energy gap at 
q = ±Q 

EM ^ S [-i?,sin7{3Ji + 23AHc) - 2J(Q)}]'/' . 

(20) 

The transverse {x, y) components of the spins arc de- 
fined in a coordinate system for which the x axis rotates 
the angle (j) around the y ot c axis from one ab plane 
to the next. In a neutron experiment with the scatter- 
ing vector k along the c axis, modulus a reciprocal lattice 
vector, the excitations detected are those described above 
at q = fc, when the scattering derives from the y ov c com- 
ponent of the spins. If the scattering is instead due to the 
X component, the experiment detects the excitations at 
q = k + Q and q ~ k — Q. The c mode behaves similarly 
to the spin waves in a simple helix, where the c compo- 
nent reflects the branch for which the Goldstone mode 
starts out from fc = 0, whereas an ab scattering vector 
component reflects the two branches, where the Gold- 
stone modes emerge from the two magnetic Bragg peaks 
at fc = ±Q. The situation is different for the w modes. 
The excitations detected by the c component are those 
determined above at g = fc, i.e. a wl 11X12 mode emerg- 
ing (subjected to a small gap) from the Bragg point at 
+Q and a w2/'wl mode emerging from the other Bragg 
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point at —Q. When the w modes are detected via the 
component in the ab plane, the wave numbers of the c- 
component branches are translated by Q or —Q. Two of 
those become pseudo Goldstone modes starting out from 
k = 0, whereas the two other branches are placed with 
their starting points at fc = 2Q and k = —2Q. The in- 
tensities of the extra ±2Q branches are, however, always 
going to be weak. 

Although the c mode shows similarities with the spin 
waves in a simple helix there is one important difference, 
namely that the cross section of the c mode vanishes when 
the total scattering vector is parallel to the c axis. In this 
case the excitations are detected exclusively via the spin 
components perpendicular to the c axis. For the c modes 
propagating along the c axis all the locally defined spin 
components within one ab plane oscillate in phase, and 
this means that the sum of the ab components for a spin- 
triangle stays zero during the c-mode oscillations. The 
two w branches starting out from ±2Q have no cross 
section either, which leaves the two pseudo Goldstone w 
branches starting out from fc = to be the only ones 
appearing in a scan along (0,0, q). Such a scan would, 
effectively, only show a single spin- wave branch (see Fig. [7] 
in the next section) , in strong contrast to that expected in 
the case of a simple helically ordered system. In general, 
dyiq) = 8y{2<j)/c-q) whereas dM ^ 3A^(l)/c-q), which 
means that there are actually two branches in the plot 
along (0,0, g) shown in Fig. [3 but the energy difference 
is small - about 0.06 meV at maximum around q — 1.25. 

In addition to the DM anisotropy included in ([3|), the 
symmetry of the system also allows the following ternii 



-Dab (fl2 • Si X S2 +f23 • S2 X S3 



fgi-SaxSi), (21) 



but this interaction docs not affect any of the quantities 
considered above to leading order in D^f^/h. The numer- 
ical analysis presented below shows that this coupling 
may weakly perturb the ordered structure by inducing a 
small oscillating c component on top of the spiraling or- 
dered moments lying in the ab plane. The ordered c com- 
ponents have a constant magnitude within an ab plane, 
which magnitude varies sinusoidally along the c axis with 
the same period as the helix. Using a D^^^ with the same 
magnitude as — 0.0038 meV considered in the next 
section, the scattering intensity due to the oscillating c- 
axis moment is found to be a factor of 10^ smaller than 
the intensity due to the a or 6 component. Hence, we may 
safely neglect any influences from D^^^. Finally, I may add 
that the simple anisotropy term D ^{S\^+ S\^ + S^^) , with 
~ 0.005 meV, may replace in the explanation for 
the confinement of the moments to the ab plane and for 
the presence of the energy gap E^{Q). In contrast to the 
DM term, the anisotropy would affect the susceptibil- 
ity components in the paramagnetic phase, but the value 
of indicated by the energy gap E^{Q) is too smaU to 
make any observable difference. However, this interac- 
tion has no influence on the choice of sign for 7, and is 
therefore unable to explain why the system prefer one of 
the two choices as found experimentallyii 
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FIG. 3: (Color online) The calculated susceptibility compo- 
nents as functions of temperature compared with the exper- 
imental results of Marty et al^ The red line and the red 
crosses denote, respectively, calculated and experimental re- 
sults, when the field is applied along the c axis. The green 
symbols denote the results obtained when the field is applied 
perpendicular to the c axis. The dashed line shows the MF 
behavior for the corresponding simple S — 5/2 system. The 
difference between the dashed and solid lines indicates the im- 
portance of correlation between the three spins in the triangle 
clusters. 



III. NUMERICAL MODEL CALCULATIONS 

In the numerical analysis, the cluster MF Hamiltonian 
is diagonalized precisely and the excitation spectra 
are calculated without making any further approxima- 
tions than the basic random-phase approximation, see 
Refs.Hl^ The best MF/RPA trimer model obtained from 
fitting the susceptibility and excitation data obtained by 
Stock et al.^ is (in units of meV): 

Ji = 1.25, J2 = 0.2, J3 = 0.1, J4 = 0.064, J5 = 0.29 
= 0.0038. (22) 

The anisotropy is very small, and I have chosen the sign 
to be positive. This means that 7 = — 27r/3, and intro- 
ducing the exchange parameter given above into the equi- 
librium condition (O leads to tancf) — 1.257 or (j> = 27r/7. 

The calculated results for the susceptibility compo- 
nents as functions of temperature are compared with ex- 
periments in Fig. [3] The dashed line in this figure shows 
the MF susceptibility for the same model given by (|22p . 
but with the intra-triangle interaction Ji included as a 
mean-field contribution. This simple MF model leads to 
a Neel temperature which is 70.1 K, whereas the present 
cluster-MF model predicts Tjy = 36.7 K in reasonable 
agreement with the experimental value = 27 K. The 
calculated Curie temperature is 6* = —190 K, and the 
one derived from the experimental paramagnetic suscep- 
tibility in Fig. |3]is 6* = —188 K (using the data within 
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Wave vector (2,0,q) (rlu) 



FIG. 4: (Color online) The calculated spin-wave dispersion 
relations along the c axis compared with the experimental re- 
sults obtained at 2.5 K by Stock et al^ along (2,0,g). The 
calculated energies of the two w modes and the two branches 
of the c mode are shown by, respectively, the blue and or- 
ange lines, and the spectrum includes only the two times two 
branches observed in the experiments. The solid lines show 
the results obtained when introducing the calculated values 
of the energy splittings and the matrix elements in Eqs. p4p 
and (|16p . whereas the dashed lines are the results obtained 
using the linear spin- wave expressions (jlSp and H19p . 

the whole interval between Tjy and 300 K), hence the 
experimental frustration factor, / = is about 7 

for this system. Because of the relatively strong intra- 
triangle interaction J^, the system is frustrated, which 
invalidates the simple MF approximation. Within the 
cluster-MF approximation, the spin triangles are treated 
as correlated units, implying that the main source for 
frustration, Ji, is accounted for in an exact way. The 
results above and the comparison in Fig. [3] show that the 
MF method is substantially improved, when choosing the 
basis to be the spin triangles instead of the individual 
spins. 

The parameters given by ([^^ lead to an exchange in- 
teraction J{Q) ~ —0.821 meV. The self-consistent di- 
agonalization of the cluster-MF Hamiltonian in ([3]) with 
h = —{S^JJ{Q) predicts the moment {gfJ-B^^z) to be 
4.71 fig/Fe in the zero temperature limit. The moment 
is reduced from its fully polarized value of 5 fj.g/Fe, be- 
cause 1 000) is not an eigenstate for the part of the 
Hamiltonian ([3]). Introducing the model parameters in 
the expression ^ for the mixing parameter, the result 
turns out to be A — 0.57. Hence, A is not small com- 
pared to 1, and the approximate expressions for the spin 
wave energies given by (dHJ and ([TO]) would not be ex- 
pected to apply. Nevertheless, the analytic spin-wave 
theory leads to useful results both when applying di- 
rectly the final results, Eqs. ^TE\\ and dT^ . or when us- 
ing instead the correct cluster MF-values for the energy 
differences and matrix elements in the expressions (|14p 



and (flBl) . as shown by respectively the dashed and the 
solid lines in Fig. 21 The numerical diagonalization of 
(1) predicts A„ = 3.486 meV and = 4.751 meV, and 
the relative matrix elements are found to be = 1.520, 
ml = 0.618, = 0.655, and = 1.275. The disper- 
sion relations along the c axis shown in Fig. 2] have been 
calculated when introducing these values in the spin- wave 
expressions dH]) and (ITB)) . The results are in good agree- 
ment with experiments and with the scattering intensity 
maxima obtained numerically. The only difference be- 
tween the numerical RPA calculations (see below) and 
the spin- wave results given by Eqs. and and 
shown by the solid lines in Fig. |6] is that the simplified 
model neglects the possible influences of the higher ly- 
ing MF levels. This minor inaccuracy is the reason why 
the calculated c-mode shows a small energy gap at Q 
and that E^(Q) is too large in comparison with the nu- 
merical RPA result. The analytical spin-wave theory is 
valuable not because it is nearly able to reproduce the 
numerical results, but because it allows a precise inter- 
pretation of the numerical RPA calculations. It is possi- 
ble to extent the analytic theory to cases where the wave 
vector also has a component in the ab plane. This is a 
more complex situation, because all three levels are being 
mixed with each other The most important change is 
that J2 in the spin-wave energies are being replaced by 
J2[cos(q • a) + cos(q • b) -I- cos(q • a-|- q • b)]/3. I have not 
found it necessary to carry through these more demand- 
ing calculations, and neither have I tried to work out 
the analytical results for the spin-wave scattering inten- 
sities, because all these extra complications are handled 
in a satisfactory and more accurate way by the numerical 
RPA calculations. 

The numerical RPA results derived for the spin-wave 
scattering intensities are presented in Figs. [5] and [51 In 
all the calculations, the intensity variations due to the 
magnetic form factor of the Fe-ions are neglected. The 
calculations of the logarithmic intensities are done with 
a narrow resolution (a Lorentzian with F ~ 0.05 meV), 
whereas a Gaussian resolution width and an intensity 
scale factor have been used as fitting parameters in the 
direct comparisons shown in Fig. [51 The two fitting pa- 
rameters are the same for all the results shown by the 
blue lines, and both parameters have been increased by 
20% in the results shown by the red lines in the right part 
of Fig. [51 These calculated results are compared with the 
experimental cross sections obtained by Stock et al^ 

As discussed in the previous section, a scan along 
(0, 0, q) should show only two nearly degenerate w 
branches starting out from the crystallographic Bragg 
point at q = 1 (in reciprocal lattice units) with the energy 
E^{Q), and the calculated scattering intensities obtained 
for such a scan are shown in Fig. [71 Preliminary measure- 
ments by Stock et al^ are consistent with this predic- 
tion. The structural parameters determined by Marty et 
al.^ show that the distance d between the ions in the 
Fe triangles is very close to be equal to (•\/3/4) a and 
thereby that b*d = n. This means that the structure 
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FIG. 5: (Color online) Logarithmic contour plots of the calculated spin-wave scattering intensities along different scattering 
vectors in reciprocal lattice units. The black squares indicate the experimental energies determined at T = 2.5 K by Stock et 
al.^ The logarithmic scale shown to the right is common for all three cases shown here and also applies to Fig. [7] and to the left 
figure in Fig. [S] The use of the logarithmic scale may be slightly misleading in the sense that the low intensity branches (light 
green/light blue), which are easily identified in these plots, are probably not visible under realistic circumstances. 
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FIG. 6: (Color online) Detail comparisons between the 
experimental^ and calculated cross sections for a selection of 
the scattering vectors considered in Fig. [5] 



factor of the spin triangles is unchanged, if 4a* or 4b* 
is added to the scattering vector. Hence, when neglect- 
ing any form factor effects, the only difference between 
a (0,0, and a (4, 0, q)/(0, 4, q) scan is that scattering 
due to the c-axis spin components, which cancels out in 
a (0,0, g) scan, contributes at (4, 0, g)/(0, 4, g). A scan 
along (4, 0, q) or (0, 4, q) is expected to show not only the 
two w-modes but also the c mode, all starting out from 
q = 1 as illustrated by Fig. [T] The w modes starting out 
from the magnetic Bragg point do (nearly) not appear 




1.1 1.2 1.3 
(0,0,q) 



1.5 1.0 1.1 1.2 1.3 1.4 
(4,0,q) or (0,4, q) 



1.5 



FIG. 7; (Color online) The logarithmic scattering intensities 
calculated along (0, 0, g) and (4, 0, g)/(0, 4, g). The experi- 
mental results included in the (4, 0, g)/(0, 4, g) figure are the 
same as those appearing in Fig. [S] and obtained by Stock et 
al^, but the results have been translated by wave vectors 

(o,o,^±g). 



in this scan, because the total c component of the spin 
triangles is zero for these modes. The experimental re- 
sults shown in the figure are the same ones as presented 
in Fig. [5l but now translated so that they all fall on the 
two branches seen in this scan. The energy gap of the w 
modes, appearing in this plot at q — 1, has been deter- 
mined in the neutron scattering experiments by Stock et 
al^ and by Loire et alA to lie between 0.35 and 0.4 meV. 
This value of the gap is here used for determining the 
numerical value of D^. This way of plotting the experi- 
mental results also shows that the energy of the c mode 
at (4, 0, l-l-Q), the upper peak in the (2,0,0) scan shown 
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in Fig. |6l is distinctively higher than the vahie suggested 
by a sinusoidal interpolation of the other c-mode results. 
The difference might be caused by special effects related 
to this particular wave vector (2,0,0), or it may be an 
indication of a weak interaction between spins in next- 
nearest neighboring layers. The fit to the two dispersion 
relations may be improved by including a ferromagnetic 
interaction Jj^jl^c) — —0.01 meV (between sublattice 1 
and 2 at a distance of 2c along the c axis), however, the 
improvements are not really significant and this possible 
modification is abandoned. 



The experiments of Stock et alA were all performed 
in the a*c* plane, whereas Loire et al.^ did choose the 
h*c* plane as scattering plane, and they observed a clear 
asymmetry between the intensities of unpolarized neu- 
trons scattered at (0,1, g) and (0, 1,— g) caused by J3 
being different from Jg. The branches starting out from 
zero energy at (0, 1,£ — Q) were found to be much more 
intense than the (0, 1,£ + Q) branches (see also the left 
figure showing 5(0, 1, q) in Fig.|5]below). This is in agree- 
ment with the theoretical predictions, if J5 is the domi- 
nant interplanar interaction, or, more accurately formu- 
lated, the theory becomes in accord with these exper- 
imental results, if the sign of R in Eq. ([7]) is chosen 
to be negative, which choice is already made with the 
model parameters introduced by Eq. (1^^ . The experi- 
ments of Loire et al^ were done on a crystal for which the 
structural chirality was determined to be e-j^ = — 1 from 
the anomalous part of the x-ray scattering function.^ 
This means that the unpolarized neutron experiments 
of Loire et al.^ show that the strongest interplanar in- 
teraction Jg is, as expectedjii^ the one determined by 
the shortest super-superexchange path. The conclusion 
that sign(i?) — ej, and therefore, according to Eq. 
EjLf = ^T^-y' or = ^H^T^ is the same one derived by 
Marty et al^ from their unpolarized neutron-diffraction 
experiments. The elastic and the inelastic unpolarized 
neutron experiments independently show that the rela- 



tion EjLf 



T "^7 



-e-P sign(Z?^) applies, but are unable 



to decide on the helicity of the magnetic structure e^j 
and, thereby, on the sign of or the orientation of the 
ordered spin triangles. The unpolarized neutron scans 
in the a*c* plane should show the similar asymmetry ex- 
cept that the high intensity branch is the one starting out 
from (1,0, Q). Preliminary results by Stock et alJ^ 
show a pronounced asymmetry between the scattering at 
(1,0,1 + Q) and (1,0,£ - Q), but the branch with the 
largest intensity is the one emerging from (1,0, — Q), 
which indicates that their crystal has = -fl, i.e. the 
opposite structural chirality to that of the crystal inves- 
tigated by Loire et al^ The coincidence that the dis- 
tance between the Fe ions in the triangles is close to 
(•\/3/4)a implies that the asymmetry disappears around 
Bragg points {h, k, t), when for example h + k = 2 like in 
the case considered in Fig. \5\ 



IV. HELICALLY POLARIZED SPIN WAVES 

Loire et alA have performed a series of inelastic scat- 
tering experiments in the h*c* plane with polarized neu- 
trons. In these experiments they determined the spin-flip 
scattering intensities I^/^ corresponding to the differ- 
ential cross sections (within the terminology of Moon, 
Riste, and Koehler^^) 



- E ^"^'^ P^Pj • Sj., T ■ (S_L, X S^^.)] 



(23) 

In their experiments Loire et al. did choose the neutron- 
spin polarization vector z to be parallel to the scattering 
vector k = kj^^j^j^^j — kg^^^^j, i.e. z = k/|k|, and extracted 
the following inelastic scattering functions^ 



5(k,L.) 



2 

/±(k,c.)-/T(k,c.) 



(24) 



In the case where the spins of the incident neutrons are 
polarized parallel to the scattering vector k, /^(k, w) is 
the intensity of the scattered neutrons with the opposite 
polarization. Since z is reversed by definition, when k is 
replaced by — k, it is found that C(k, uj) = C(— k, w) as 
well as ^(k, uA — iS'(— k, uj). The expression given by Eq. 
(16) in Ref. [a is straightforwardly generalized so to in- 
clude C(k, cj): the products of spin components appear- 
ing above in Eq. p3)) are translated into the correspond- 
ing tensor components of the correlation function in Eq. 
(16) of Ref. 0, but notice that the sign convention for the 
wave vectors applied in this reference is the opposite of 
that used above. For a simple righthanded (e^ = +1) 
helix with ordering wave vector Qc* (0 < Q < 1/2), 
the elastic response P is non zero, and is zero, for 
the magnetic Bragg peaks at k = (0, 0, \£\ + Q) and 
(0, 0, —\£\ — Q), whereas the opposite is the case for the 
remaining Bragg points along c* . The ratio between the 
static scattering functions C(k) and ^(k) for a simple 
helix with helicity is in general found to be 



2 cos 6* 



g(k) 

S'(k) 1 + cos2 ( 



[5{G - Q - k) - <5(G + Q - k)] e^, 

(25) 

where cos6' = k • Q/|k||Q| and G is a reciprocal lattice 
vector. 

The polarized inelastic neutron-scattering experiments 
performed by Loire et alf^ on a crystal with e^^ = — 1 
showed that C(k,CLi) is positive and very nearly equal 
to ^(k, a;) at the lowest energies observed (1 meV) close 
to the Bragg point at (0,1,6/7), and that this is also 
the case for the c mode at (0,1,1.3) as shown in Fig. 
m The behavior of C(k, w) was measured along (0, 1,<?), 
and the results shown by Loire et al. in their Fig. 3(e) 
are very similar to the calculated results shown in the 
middle figure in Fig. [5] (the weak branches included here 
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FIG. 8: (Color online) Contour plots of the calculated spin-wave scattering intensities along k = (0, !,<?). The left figure 
shows the logarithmic intensities S(k, lj) derived for unpolarized neutron scattering. Notice that the intensities of the branches 
starting out from zero energy at (0, 1, £ — Q) are much higher than the intensities of the {0,l,l + Q) branches. The figure in the 
middle shows the calculated polarization function C(k, lu) using a linear scale. The figure to the right shows the corresponding 
dynamic scattering ratio C(k,ijj)/ S(k,uj). 
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FIG. 9; (Color online) Comparison between theory and ex- 
periments for two of the scattering vectors considered in Fig. 
[8] The experimental results are obtained by Loire et al,"^ from 
inelastic scattering of polarized neutrons. The intensity scal- 
ing factor and the resolution width are kept constant in all 
the calculated results. 



were not detectable). The sign of C(k, w) determines 
the sign of the helicity of the helically ordered moments, 
and, as concluded by Loire et al.f" their results show 
that the helicity — +1. This result, in combination 
with the unpolarized result ej, — sign(i?) = — 1 and the 
maximization of 3(Q), Eq. ([8]), implies that = —1 
corresponding to a positive value of D^. 

For a simple helix, the spin-wave theorj*i^ii^ pre- 
dicts the dynamic ratio between the scattering functions 
C(k, a;)/S'(k, w) to be a purely geometric quantity, like 
the static one in Eq. (PS)) . 



C(k,w) _ 2 cos 6* 
S{k,uj) ^ ^ 1 + cos2 e 



(26) 



This expression applies to the modes emerging from the 
magnetic Bragg points, which scatter the neutrons via 
the planar spin components. The sign to be chosen for 
the ± sign in front is the same as the sign in front of 
the corresponding delta function in the magnetic Bragg 
scattering ratio in Eq. (P5|) . i.e. the -I- sign applies to 



the modes starting out from the Bragg points at G — Q. 
In contrast to this, the cross section deriving from the 
c component of the spins are not affected by the helic- 
ity, and C(k, cij) — for the branches emerging from the 
nuclear Bragg peaks. With the restriction that we shall 
only consider spin waves propagating along the c axis, 
this simple description applies almost unchanged to the 
present ordered system consisting of three sublattices of 
coupled helices. C(k, w) = for the c modes emerg- 
ing from the nuclear Bragg peaks at G and for the w 
modes starting out from the magnetic Bragg peaks at 
G ± Q. The c modes detected via the spin components 
in the ab plane, are the branches emerging from the mag- 
netic Bragg points G ± Q, and the chiral scattering ratio 
for these modes is determined by Eq. (|26p . In principle, 
the w modes starting out from a nuclear Bragg point 
or from G ± 2Q should show a dynamic scattering ra- 
tio which is also determined by Eq. ((26)) . However, the 
±2Q branches are weak, and the two w modes emerging 
from the same Bragg point are not easy to separate, and 
since the chiral scattering ratios have opposite sign for 
the two branches, they are going to appear like a single 
mode with a chiral polarization factor close to zero. In 
the scan along (0, 0, q) shown in Fig. [71 the upper one of 
the two nearly degenerate w modes, the —Q branch, has 
C{k,uj)/S{k,u}) = +1 whereas C(k,Lj)/5(k,Lj) = -1 for 
the other +Q mode. These results are valid in the case 
where the spin waves are propagating along the c axis, 
and they agree in most details with the dynamic scatter- 
ing ratios calculated numerically using the RPA model. 
One example is shown in the right of Fig. [51 In prin- 
ciple, the colors in this figure should become more and 
more blue [C (k, uj) / S (k, uj) — +1] or more and more red 
[C{k,u!)/ S{k,u!) = —1] for increasing values of q, how- 
ever, this systematic behavior is going to be disturbed 
whenever intensities from modes with different helicity 
factors overlap each other (notice that \C{k,uj)/ S{k,uj)\ 
is predicted to be 0.967, nearly 1, already at q = 6/7). 
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The inelastic cross section is readily calculated numer- 
ically in the general case, but the result becomes less 
transparent due to the complication that the c and w 
modes are coupled whenever the propagation vector has 
a non-zero component in the ab plane. 

V. CONCLUSION 

The intra-triangular interaction Jj^ > is the domi- 
nant cause for the relatively strong frustration shown by 
the present spin system, and it is important to account 
for this interaction primary to the interactions between 
the Fe triangles. The MF approximation, when applied 
to the trimerized clusters, leads to a much improved de- 
scription of the system in comparison with the single-spin 
approximation. The susceptibility just above is be- 
ing reduced by 15% and the transition temperature by 
50%. The more precise description of the ground state 
is essential for the determination of the exchange con- 
stants. The analysis of the spin waves in terms of the 
five exchange constants, plus the DM anisotropy, leaves 
one exchange constant as a nearly free variable. In the 
present analysis this degree of freedom has been fixed 
by a fitting to the susceptibility. The present model ac- 
counts very well for the susceptibility, whereas the use 
of, for instance, the exchange parameters proposed by 
Loire et al^ leads to values for the susceptibility, which 
are about 25% larger than the experimental ones (40% in 
the single-spin MF approximation) at temperatures close 
to T^. The MF-cluster calculation also predicts a quan- 
tum reduction of the ordered moment, from 5 to 4.7 /i^. 
This reduction is less pronounced than the one derived 
from experiments, as the value of the ordered moment 
determined from neutron-diffraction experiments at 2 K 
by Marty et alX was found to be as small as about 4 /i^ . 

Although the trimerized states of the Fe spin trian- 
gles differ greatly from that determined by fully polarized 
spins, the linear spin- wave theory works surprisingly well 
demonstrating the effectiveness of the boson represen- 
tation obtained by, for instance, the Holstein-Primakoff 
transformation. This advantageous property is of similar 
significance for the analysis of an anisotropic ferromag- 
nct like terbium metal, which is discussed at length in 
Refs.lsllTsl. The particular model derived by Loire et al^ 
on the basis of linear spin- wave theory is questioned, but 
their theoretical results are found to agree closely with 
the predictions of the present theory. Their characteriza- 
tion of the different spin-wave modes may be imprecise, 
but their main conclusions are the same as derived here. 
In the crystal they investigated the crystallographic chi- 
rality is = —1, the ordered moments in the triangles 
have the orientation specified by e..^, ~ — 1 , and the spiral- 
ing ordered moments along a line parallel to the c axis are 



making a righthanded helix corresponding to = -1-1. 

The mirroring of the present system with respect to, 
for instance, the ac plane through Sj^ would change the 
structural chirality e-^ from —1 to -1-1, and the helicity 
e^j from -|-1 to —1, whereas the sign of D^{> 0) and 

= —1 would remain unchanged. The choice between 
the two possible orientations of the ordered spin trian- 
gles is determined by the local surroundings not by the 
structural chirality. By definition chirality should be in- 
variant with respect to time reversal but have odd parity 
with respect to inversion. Based on this definition, the 
orientation of the spin triangles denoted by cannot be 
characterized as a chiral property. 

To summarize: The asymmetry between the inelastic 
intensities of magnetic scattered unpolarized neutrons at 
(1, 0, i-q) and (0, 1, l+q) or at (1, 0, l-q) and (1, 0, l+q), 
as observed, respectively, by Loire et al^ and by Stock 
et ali^, is alone a consequence of the structural chirality. 
ct is equal to —1 in the crystal studied by Loire et ai, 
whereas e^^ = -|-1 is the preliminary result for the crystal 
investigated by Stock et al. The DM anisotropy term 
in Eq. ([2]), due to the spin-orbit interaction, couples 
chirality in spin space to the crystallographic chirality, 
ejj = e^erp. The sign of the DM anisotropy determines 
the choice between the two possible orientations of the 
ordered moments in the Fe triangles, and it is found that 
= — sign(Z?^) = —1. The combination of an enan- 
tiopure crystal and a non-zero DM anisotropy implies 
that the system only contains ordered moments with a 
single sense of helicity, that = —ej, in this system 
where is positive. The dynamics of the present sys- 
tem is unique, not because the spin waves have chiral 
properties but because the presence of only one domain 
of helicity has made it possible to observed this intrinsic 
dynamic property. It may be concluded that the three 
different modes of spin waves propagating along the c axis 
in the Ba3NbFe3Si20i4 crystal investigated by Loire et 
alA should all posses, depending on their effective propa- 
gation vector q = k — G, the same or the opposite sense 
of helicity as the ordered structure. The observed behav- 
ior of the c mode is in agreement with this conclusion, 
whereas the near degeneracy of the w modes starting out 
from the nuclear Bragg points prevents an experimental 
determination of the chiral properties of these modes. 
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